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Abstract 



The reactions of cold H atoms with solid O2 molecules were investigated at 10 K. The 
formation of H2O2 and H2O has been confirmed by in- situ infrared spectroscopy. 
We found that the reaction proceeds very efficiently and obtained the effective 
reaction rates. This is the first clear experimental evidence of the formation of water 
molecules under conditions mimicking those found in cold interstellar molecular 
clouds. Based on the experimental results, we discuss the reaction mechanism and 
', astrophysical implications. 
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1 Introduction 



Water is the most abundant solid material in space, and has been observed in 
various astrophysical environments, such as outer planets, satellites, comets, 
and interstellar clouds pQ. Since the solar system evolved from an interstellar 
molecular cloud, icy objects in the solar system originated from the water ice 
formed in the interstellar molecular cloud. Therefore, gaining an understand- 
ing of the origin of water molecules in interstellar molecular clouds is critical 
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not only for discussing the origin of the solar system, but also for under- 
standing chemical evolution and the origin of life pj. However, the formation 
mechanism of water molecules in the interstellar clouds has not been under- 
stood to date. It has been clarified that the formation of water molecules in 
the gas phase is incapable of explaining the observed abundance in molecular 
clouds [3114"] . Thus, it has been suggested that water molecules are synthesized 
by atomic reactions involving H and O on pre-existing silicate or carbonaceous 
grains at around 10 K [31I4T5] . 
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Fig. 1. Possible reaction scheme for formation of H2O involving H and O atoms in 
interstellar clouds (modified from |5|8j ). Atoms or molecules and numerals next to 
the arrows indicate the reactants in the reaction and the number of the reaction in 
the text, respectively. Broken and solid arrows indicate reactions with and without 
activation barriers, respectively. Thick and thin arrows denote reactions investigated 
and not investigated in the present study, respectively. Molecules observed in the 
present study are enclosed in thick squares. 



The first possible process for H 2 formation is the sequential hydrogenation 
of O atoms on grains: 

O + H^OH, (1) 

OH + H^H 2 0. (2) 

These reactions have no activation barriers j6] The second possible process is 
the reaction of OH with hydrogen molecules absorbed on the surface of grains: 



OH + H 2 -> H 2 + H, 



(3) 



which has an activation energy of 5.16 kcal/mol in the gas phase[7]. For the 
third process, Tielens and HagenjSj proposed alternative reactions: 



O + O^O 



2- 



2 + H -> H0 2 , 
H0 2 + H^H 2 2 , 
H 2 2 + H -> H 2 + OH. 



(4) 

(5) 
(6) 
(7) 



Reaction (jSJ) has essentially no barrier; theoretically estimated activation ener- 
gies lie between 0.1 and 0.4 kcal/mol [9]. Reaction (JTj) has activation energies of 
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3.6-4.3 kcal/mol [TU] . All activation energies referred here are value in the gas 
phase. There has been no data available on a surface. Fig. [1] summarizes the 
possible water formation routes from reactions (fTT)-(r7j). Cuppen and Herbst [5] 
simulated the formation of water molecules in various environments of inter- 
stellar clouds, and found that reactions ([1]) and ([2]) are the main routes in cool 
diffuse clouds and reactions from ([3])- ([7]) are the main routes in cold molecular 
clouds. H and O atoms are the major gas-phase species in the former clouds, 
whereas H 2 molecules are the major gas-phase species in the latter clouds. 

Hiraoka et al. [TT] investigated reactions ([T]) and §Z§ by spraying D atoms onto 
O atoms trapped in an N 2 matrix at 12 K. They analyzed the products by 
using temperature-programmed desorption (TPD) spectroscopy and observed 
D2O. A similar experiment has been performed by Dulieu et al. [12] although 
their results are still very preliminary. However, in the experiments of both 
groups, it remains unclear whether D 2 is formed at 10 - 12 K or during 
heating. It is essential to perform in- situ analysis during H/D atom irradiation 
to confirm the formation temperature of H 2 and to obtain kinetic data. For 
reactions ([S])-©, Klein and Scheer [13] performed pioneering experiments on 
the reaction of H atoms with solid 2 at 20 K. Although they speculated on 
the production of H 2 2 and H 2 based on their experimental results, they 
were unable to verify the production of these molecules since no analysis of 
the products was performed. To investigate the formation mechanism of water 
molecules in interstellar molecular clouds, we focus on reactions ©-([7]) and the 
subsequent reactions ([2]) and ([3]), and performed H/D addition experiments 
to solid 2 at 10 K. 



2 Experimental 

Experiments were performed using the Apparatus for SUrface Reaction in 
Astrophysics (ASURA) system described previously [14715] . Briefly, solid 2 
with an 8 monolayer (ML) thickness was produced by vapor deposition on an 
aluminum substrate at 10 K at a deposition rate of 25 pm s -1 in an ultrahigh 
vacuum chamber at a pressure of about 10~ 8 Pa. Atomic H and D were pro- 
duced by a microwave-induced plasma in a Pyrex tube and transferred via a 
series of PTFE and aluminum tubes to the target. The dissociation fraction 
was about 20% at least. The atomic beams could be cooled to 20 K in the 
aluminum tube that was connected to an He refrigerator. The flux of atoms 
was controlled by varying the temperature of the aluminum tube and mea- 
sured by the method of Hidaka et al. [16J. The fluxes of both the H and D 
atoms were 2xl0 14 cm -2 s -1 in the present experiments. Infrared absorption 
spectra of the sample solid during irradiation by atoms were measured in-situ 
by Fourier transform infrared spectroscopy with a resolution of 4 cm -1 . After 
irradiation, TPD spectra of the sample were obtained using a quadrupole mass 
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spectrometer at a heating rate of 4 K min 



3 Results and discussion 
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Fig. 2. Spectral change in solid O2 due to (a) cold H exposure and (b) D exposure. 
The spectra were obtained by subtracting the initial absorption spectra of non-ir- 
radiated samples from the spectra of H(D)-irradiated samples. Lines indicate newly 
formed molecules. Spikes with arrows indicate noise caused by vibration of the He 
refrigerator. 



Fig. 2a shows the change in a typical absorption spectra on irradiation by 70- 
K-hydrogen atoms at 10 K. This figure shows the absorbance variations from 
the initial spectrum of the sample, so that bands appearing above the baseline 
indicate an increase in absorbance. H2O2 appears immediately with bands at 
3250, 2830, and 1405 cm" 1 , then H 2 follows at 3432, 1650, and 820 cm" 1 . 
No product was observed when solid O2 was exposed to an H 2 beam at 10 



(a) 



t = 5 sec 




\ Noise 




K, or when H atoms were sprayed onto the Al substrate without solid O2, 
demonstrating that the formation of H2O2 and H 2 in the present experiments 
is due to H atom irradiation onto O2. No intermediate radicals (e.g., HO2) were 
observed, suggesting that the reaction rate of ([SD is much slower than that of 
([6]). In this experiment, the H flux is sufficiently large so that the reaction 
(jSJ) immediately proceeds to produce HO2. The formation of H2O2 and H 2 
were also observed when solid O2 was exposed to an H beam at temperatures 
between 15 K and 20 K. Fig. 2b clearly reveals the production of D2O2 and 
D 2 in D-exposure experiments, confirming that the H 2 produced in the 
present experiments is not a contaminant. The production of H 2 02(D 2 02) and 
H20(D20) was also indicated in TPD spectra obtained after H(D)-exposure. 

The amount of products are calculated and plotted as a function of expo- 
sure time in Fig. 3. These were obtained from the band areas in the spectra 
and the reported and estimated integrated band strengths. The integrated 
band strengths used are 1.2xl0~ 17 and 2.1xl0~ 17 molecule cm -1 for H2O 
(1655 cm -1 ) [17J and H2O2 (1347 cm -1 ) [18], respectively. Since no integrated 
band strengths for D 2 at 1200 cm -1 and D2O2 at 1033 cm -1 have been 
reported, we estimated these in the following manner. Bergren et al. [19] mea- 
sured the optical constants of OH(OD)-stretching modes in amorphous H 2 
and D2O at 3200 cm -1 and 2400 cm -1 , respectively. Assuming that the ratio of 
integrated band strengths for H2O and D2O in OH(OD)-stretching modes are 
the same as those in OH(OD)-bending modes, we estimated the integrated 
band strength of OD-bending modes for D 2 (1200 cm -1 ) to be 8xl0 -18 
molecule cm -1 . No data is available for D2O2, not even optical constants. So 
we assumed that the ratio of integrated band strengths for H20/D20= 0.67 is 
the same as that for H2O2/D2O2, and estimated the integrated band strength 
of OD-bending modes for D2O2 (1033 cm^ 1 ) to be 1.5xl0~ 17 molecule cm -1 . 
This result is consistent with the spectral features of OH and OD-bending 
modes of amorphous H2O2 and D2O2 measured by Lannon et al. [20], which 
show that the ratio of peak heights for H2O2/D2O2 is 0.75. 

Fig. 3a clearly shows that the formation of H2O2 and H2O is very rapid and 
efficient; H2O2 and H 2 are observed even after exposure for 5 s. When H 
atoms are irradiated onto solid CO at 10 K [T4], the production of H 2 CO and 
CH3OH are observed after 30 s and 1 min, respectively. The present exper- 
imental results confirm that reactions which were initially proposed 
based on theoretical considerations 06], proceed at 10 K. Since reaction (jSJ) 
essentially has no barrier [9J, it is natural that this reaction proceeds very 
rapidly. While reaction (j7j) has an activation energy of 3.6-4.3 kcal/mol JTU], it 
has been suggested that this reaction does not proceed by an Arrhenius-type 
reaction at 10 K, but rather proceeds by a tunneling reaction even at 10 K. 
As described later, isotope effect (H or D) on reaction ([7]) has been observed, 
strongly supporting that reaction ((Tj) proceeds by the tunneling reaction. Af- 
ter the formation of OH by reaction fl7|), there are two possible reaction routes 
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Fig. 3. Variation in column densities of products (a) for H2O2 and H2O with H 
exposure and (b) for D2O2 and D2O with D exposure. Each data point represents 
the average of three measurements, and error bars represent the statistical errors. 
The upper abscissa indicates the fiuence of H(D) atoms. The solid lines are fitted to 
the curves defined by Eqs. (fTT|) to (fT3|) . The amount of O2 reduction is assumed to 
be equal to the sum of the increase in H2O2 (D2O2) and I/2H2O (I/2D2O). Fitting 
curves of H2O and D2O shown are two times larger than those obtained from Eq. 

for further H2O formation: reactions (j2j) and ([3]). This is because the H-beams 
used in the present experiments do not consist of pure H atoms, but rather 
are a mixture of H and H2. Reaction (j2J) has no activation barrier, while the 
activation energy of reaction ([3]) is 5.16 kcal/mol [7] for reaction ([3]). It seems 
probable that H 2 is formed by reaction (j2j), although we currently have no 
direct evidence for this. Fig. 3b shows that D 2 2 and D 2 are produced very 
rapidly and efficiently in the D-exposure experiments, although the formation 
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of D2O is slightly slower than in H-exposure experiments. 



In the present experiments, it is impossible to directly measure the decrease 
in O2, since O2 is infrared inactive. Thus, it is reasonable to assume that the 
summation of 0.5H 2 O and H 2 2 is equal to the reduction of 2 . The 2 loss 
becomes saturated after long exposure. This type of saturation was also found 
previously for the successive hydrogenation of CO at 8-15 K by Watanabe et 
al. [Tip!] . They concluded from the dependence of the saturation value on the 
sample thickness that the saturation is related not to the balance between the 
forward and reverse reactions, but rather to the diffusion length of H in solid 
CO and CO-H 2 mixed ice. This suggests that in the present experiments 
the reaction occurs only at the surface. The structure of the molecular layer 
may thus be onion-like, with H 2 on the top, H 2 2 in the middle, and 2 
at the bottom. The amount of H 2 2 is greater than that of D 2 2 , suggesting 
that the diffusion length of H into 2 is longer than that of D into 2 . On 
the other hand, the amount of H 2 is the same as that of D 2 0, suggesting 
that the water formation reaction occurs only at surface. The small amount of 
observed H 2 compared to H 2 2 may be attributed to desorption of products 
caused by the release of the heat of reaction. 

To discuss the reaction mechanism based on the isotope effect, we next derive 
effective reaction rate constants by fitting the data in Fig. 3. To simplify the 
problem, we assume that reactions ©-([Tj), ([2]) and (j3j) can be simplified to 
the following consecutive reaction: 



where k{ is the reaction rate constant of the ith reaction. Since, it is very 
difficult to measure the surface density of H atoms, nn, in the present study, we 
were not able to obtain ki directly. Instead, we define effective rate constants : 
k' m = kiUii. Since n H is considerd to be time- independent, these reactions are 
assumed to be pseudo-first-order reactions. 

Furthermore, we assume that the products have a layered structure, as men- 
tioned above. In the deepest layer, no reaction of 2 with H atoms occurs. 
In the middle layer that has a thickness of a, the reaction from 2 to H 2 2 
occurs, but that from H 2 2 to H 2 does not occur. In the uppermost layer 
with a thickness of b, the entire reaction (|S]) occurs. The rate equations for the 
number densities of 2 , H 2 2 and H 2 are given by 



2 



H 2 2 



2H 2 
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The solutions of equations Q)- (TTTT) for the initial conditions of n (t = 0) = n, 
and ni(0) = n 2 (0) = are given by 

n 



be fcm * + ae k ' ni 



n° 



m _ 








"-Hi 


- k 


n 2 /2 
"o 


= b 


1 + 



y y \ m H1 

^Hl "-H2 



(12) 

(13) 
(14) 



where Uq + rii + n 2 /2 = rip. In the case of D exposure, the same reaction as 
that for H exposure is assumed to occur, namely, 



2 ^ D 2 2 D 2 0, 



(15) 



where fc^ = A^riD, and ra^ is the surface density of D atoms. Effective rate 
constants, k' m and k' Di , obtained by fitting the data are shown in Table 1 
together with those for the first steps of successive hydrogenation (deuteration) 
of CO [15]fT6] : 

CO ^ HCO -> H 2 CO -» CH 3 -> CH3OH, (16) 

CO ^ DCO -> D 2 CO -> DH3O -> DH3OD. (17) 

Using the relation of %/nD = 1 [16] . the ratios of the reaction rate constants 
are estimated: k' m /k' m = 1 and k' m /k' D2 = 8. The former result states that 
there is no difference in rate constants between 2 + H and 2 + D at 10 K, 
and it is consistent with the fact that reaction (JSJ) has essentially no barrier [5]. 
The latter result for the ratio of reaction rate constants between H 2 2 + H 
and D 2 2 + D {k'^/^m = 8) is reasonable for a tunneling reaction with an 
activation energy of 3.6-4.3 kcal/mol flOj. In the case of CO + H(D), the 
ratio of reaction rate constants, k^/k'^, lies between 10 (from Table 1) and 
13 [TB] with an activation energy of about 4 kcal/mol. In any case, these data 
are very useful for estimating the barrier height and width of a potential, 
and potential energy surface as discussed by Hidaka et al. [16J. As already 
mentioned qualitatively, the rate constants of (jHJ) and (fT5|) are one to two 
orders of magnitude greater than those of (fTB]) and (TTTT) . 



4 Astrophysical implications 



Our results have several implications for ice in space. Assuming a hydrogen 
number density of 1 cm -3 , the H fluences in a 10-K molecular cloud will be 
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Table 1 



Effective reaction rate constants (min 1 ) 



Temperature 






Reaction 








(K) 


2 +H 


H 2 2 +H 


2 +D 


D 2 2 +D 


CO+H 


CO+D 




k' 


k' 


k' 


k' 


k' 


k' 


10 


12.8 


3.9 


12.0 


0.49 


0.14 a 


0.014 fe 


15 










0.41 c 


0.033 c 



a,b,c: obtained with the H(D) flux of lx 10 14 cm- 2 s -1 [22 



a: for pure CO |15j : b: for pure CO |16j ; c: for 0.8 ML CO on amorphous H 2 [16] 

1.3 xlO 16 and 1.3 xlO 17 cm -2 over 10 4 and 10 5 years, respectively. In the present 
experiment, these fluences approximately correspond to exposure times of 1 
and 10 minutes. Considering that the lifetime of molecular clouds is of the 
order of 10 6 years, it is reasonable to say that water formation by reactions (jSJ)- 
(J7|) occurs very quickly in molecular clouds. Based on computer simulations of 
surface reaction networks, Tielens and Hagen [8] and Cuppen and Herbst (6] 
conjectured that reactions (J1])-(J7]) occur in cold molecular clouds in which 
the ratio H/H 2 is low. The present results strongly support their conclusion 
experimentally, although the contribution of reactions ([I])-© is still unclear. 

Solid 2 is expected to occur on the grain surface with H 2 and CO ice in 
interstellar molecular clouds and to be observable by infrared telescopes since 
the 6.44-/xm 0=0 stretching mode is observable when 2 is surrounded by 
other molecules [23J and because the 4.67-//m feature of CO is altered when 
2 is mixed with CO [21]. However, no positive detection has been reported 
so far [21] . Present experiments clearly explain why no 2 has been observed: 
if an 2 molecule is formed on grain surface it will react very quickly with H 
to form H 2 2 and H 2 within 10 4-5 years as shown in Figs. 3a. 

Recent astronomical observations have revealed that the abundances of some 
deuterated molecules are up to four orders of magnitude greater than the cos- 
mic D/H ratio of 1.5xlCT 5 . In the case of water, a HDO/H 2 ratio of 0.03 
has been measured in a solar-type protostar [25j. Many theoretical models, 
including pure gas phase models [2S] and gas-grain models [27], have been 
proposed. However, only a few experimental studies have been performed to 
verify the latter models [16f28] . In the case of formaldehyde and methanol, 
the successive addition of D to CO (reaction fflTj) ) is not favorable for produc- 
ing deuterated formaldehyde and methanol [IB] . Instead, H-D substitution in 
formaldehyde and methanol are necessary to achieve the observed abundances 
of deuterated spiecies [2E] ■ However, Nagaoka et al. [28] found experimentally 
that no deuteration of H 2 occurs by H-D substitution under D exposure at 
10-20 K, even for fluences of up to 5xl0 18 cm -2 . Considering the ratios of 
k' m /k' D1 = 1 and kj 32 /k^ 2 = 8 obtained in the present experiments and the 
D/H atom ratio of 0.1 or less expected in molecular clouds [2H], deuterium 
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addition to O2 is favorable for producing the observed amount of HDO; the 
observed HDO/H 2 ratio of 0.03 can be achieved in a time scale between 
10 4 and 10 5 years. Although the discussion is indefinite at present, the ex- 
perimental results presented in this study provide a basis for discussing H/D 
fractionation. 
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